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We report numerical evidence that an epidemic-like model, 
which can be interpreted as the propagation of a rumor, ex- 
hibits critical behavior at a finite randomness of the under- 
lying small-world network. The transition occurs between a 
regime where the rumor "dies" in a small neighborhood of its 
origin, and a regime where it spreads over a finite fraction 
of the whole population. Critical exponents are evaluated, 
and the dependence of the critical randomness with the net- 
work connectivity is studied. The behavior of this system as 
a function of the small-network randomness bears noticeable 
similarities with an epidemiological model reported recently 
[M. Kuperman and G. Abramson, Phys. Rev. Lett. 86, 2909 
(2001)], in spite of substantial differences in the respective 
dynamical rules. 

PACS numbers: 87.23. Ge, 89.75.Hc, 05.10.-a 



Small-world networks have been introduced as an in- 
terpolation between ordered and random graphs, to cap- 
ture two specific features of real neural, social, and tech- 
nological networks Jl|,||. On one hand, they have a rel- 
atively large clustering coefficient, i.e. the probability 
that two neighbors of a given vertex are in turn mu- 
tual neighbors is high, as in ordered graphs. On the 
other, the typical separation between two vertices is much 
smaller than the total number of vertices, as in random 
graphs. Whereas the geometrical properties of small- 
world networks have been studied in detail £0-9], less 
attention has been paid to the dynamical properties re- 
sulting from their partially random structure P|MLll] . As 
explained in more detail below, the structure of small- 
world networks is parametrized by the randomness p, 
where p = and p = 1 correspond to fully ordered and 
random neworks, respectively. It has been shown that, 
in asymptotically large small- world networks, statistical 
geometrical properties — such as the mean distance be- 
tween vertices — exhibit qualitatively the same behavior 
for any p > 0. Here, we study an epidemic-like propaga- 
tion process on a small-world network which, in contrast 
with geometrical properties, shows a transition between 
two qualitatively different dynamical regimes at a finite 
value of p. Since this model addresses a social process 
and small-world networks are a presumably realistic rep- 
resentation of social networks, the analysis could be rel- 
evant to the description of threshold phenomena in real 
societies. 

Our model [ ^2|JT^ ] consists of an A^-element popula- 
tion where, at each time step, each element adopts one 



of three possible states. By analogy with epidemiological 
SIR models these states are called susceptible (S), 
infected (I), and refractory (R). The evolution proceeds 
as follows. At each time step a randomly chosen infected 
element i contacts another element j. Then, (i) if j is in 
the susceptible state, it becomes infected; (ii) if, on the 
contrary, j is infected or refractory, i becomes refractory. 
These rules are better interpreted in the frame of a rumor 
spreading process, where S-elements have not heard the 
rumor yet, I-elements have heard the rumor and are will- 
ing to transmit it, and R-elements have lost their interest 
in the rumor and do not transmit it. Initially, only one 
element is infected and the remaining N — 1 elements are 
susceptible. During the first stage of the evolution, the 
number of I-elements increases. Since this also implies 
a growth of the R-population, the contacts of I-elements 
between themselves and with R-elements become more 
frequent. After a while, in consequence, the I-population 
begins to decline. Eventually, it vanishes and the evo- 
lution stops. At the end, Nr elements — now in the re- 
fractory state — have been infected at some stage during 
the evolution. Numerical simulations and analytical re- 
sults show that, generally, Nr < N. Therefore, there 
is a fraction of the population that never hears the ru- 
mor. In the original version of this model, each I-element 
is allowed to contact at random any other element of 
the population. It has been proven that, in such situ- 
ation, the ratio Nr/N approaches a well-defined limit, 
Nr/N — 0.796 . . ., for asymptotically large values of N 

IIS 

Here, in contrast, we assume that the elements are 
situated at the vertices of a small-world network, and 
contacts can be established between linked elements 
only. The small-world network is constructed from a 
one-dimensional ordered network with periodic bound- 
ary conditions — a ring — where each node is linked to its 
2K nearest neighbors, i.e. to the nearest K neighbors 
clockwise and counterclockwise p|p|Jri|]. Then, each of 
the K clockwise connections of each node i is rewired 
with probability p to a randomly chosen node j, not be- 
longing to the neighborhood of i. A short-cut between 
two otherwise distant regions is thus created. Double and 
multiple links are forbidden, and realizations where the 
small- world network becomes disconnected are discarded. 
As advanced above, the parameter p measures the ran- 
domness of the resulting small-world network. Note that, 
independently of the value of p, the average number of 
links per site is always 2K. We have performed series of 
10 3 to 10 5 numerical realizations of the model for several 
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values of p, N, and K. At each realization, the small- 
world network was generated anew and the evolution was 
recorded until the exhaustion of the I-population. 
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FIG. 1. Frequency distribution of the number Nr of 
R-elements at the end of the evolution, for K = 2 and (a) 
p = 0.05, (b) p = 0.19, and (c) p = 0.3. Different symbols 
correspond to N = 10 3 (crosses), N = 10 4 (empty dots), and 
N — 10 5 (full dots). In (b), the dotted straight line corre- 
sponds to a power law ~ Nr ■ The dotted curve in (c) is 
a schematic representation of the iV-independent profile ob- 
served for small Nr. Note carefully the different scales of the 
three plots. 

In the first place, we have studied the distribution of 
the number Nr of R-elements at the end of the evolution 
for K = 2. Figure [l] shows the frequency /(Nr) ob- 
tained from series of 10 5 realizations for selected values 
of p and N. For small randomness (p = 0.05, Fig. |l]a) 
the distribution is approximately exponential and does 
not depend on N. In this regime, the rumor "dies" after 
a few time steps in a small neighborhood of its origin, 
due to the lost of interest of the highly interconnected el- 
ements close to the initially infected site. Therefore, the 
size N of the whole population is irrelevant to the value of 
Nr. The situation is considerably different for relatively 
large randomness (p = 0.3, Fig. |l|c). Here the distribu- 
tion /(Nr) is bimodal, with a maximum close to Nr = 



(not shown in the figure) and an additional bump for 
larger Nr. Near Nr — 0, the frequency is independent 
of N, as in the case of small randomness. Contributions 
to this zone of the distribution come from realizations 
where propagation ceases before a short-cut is reached. 
In contrast, the additional structure is strongly depen- 
dent on N. The position of its maximum, in fact, grows 
linearly, as iVp lax w 0.25./V. In a typical realization con- 
tributing to this zone of the distribution, many contacts 
occur through short-cuts and a finite portion of the pop- 
ulation becomes infected. The intermediate regime, just 
before the large-A^ structure begins to build up, is il- 
lustrated in Fig. |l|b for p = 0.19. The frequency follows 
here a power law, /(Nr) ~ N^ a with a w 1.5, over a 
substantial interval. This interval is limited by above by 
a smooth cut-off, which shifts to larger values of Nr as 
N p , with P ~ 0.5. 
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FIG. 2. Average fraction r — (Nr/N) of refractory ele- 
ments at the end of the evolution as a function of the ran- 
domness p on a small- world network with K — 2, and for 
N = 10 3 (dots), N = 10 4 (dashed line), N = 10 5 (dotted 
line), and N = 10 6 (full line). The insert shows a close-up in 
the transition zone. The arrow indicates the critical value of 
p, determined as explained in the text. 

The appearance of a well-defined power-law distribu- 
tion is a clue to the critical-phenomenon nature of the 
transition from the regime where the rumor remains lo- 
calized to the regime where it spreads over a finite frac- 
tion of the population. To characterize the transition, we 
choose as an order parameter the average fraction Nr/N 
of R-elements at the end of the process, 

N 

r = (N R /N) = N- 1 J2 N R f(N R ), (1) 

N R =Q 

which can be straightforwardly calculated from the nu- 
merical results. In the small-randomness regime, where 
J(Nr) is independent of N, we expect r ~ iV -1 , so that 
r — > as N — > oo. For large p, in contrast, the pres- 
ence of the large- Nr maximum in /(Nr) should provide 
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a finite contribution to r even for asymptotically large 
populations. These features are verified in numerical re- 
alizations. Figure || shows r as a function of p, calculated 
in series of 10 3 to 10 4 realizations for N ranging from 10 3 
to 10 6 and K = 2. A well-defined transition at a finite 
randomness p c is apparent for large N. Assuming criti- 
cal behavior at the transition, p c can be determined as 
the randomness for which the extension of the power- 
law interval in /(TVr) is maximal. This criterion yields 
p c = 0.19 ±0.01. 

The insert of Fig. ^ shows a close-up of the main plot 
near p c . The smoothness of the curves, observed even 
for the largest values of N, suggest that the critical ex- 
ponent 7 associated with r ~ \p — p c \ 7 just above the 
transition, should be larger than unity. Linear fitting of 
r as a function of p — p c in a log-log representation, for 
N = 10 5 and 10 6 , gives 7 » 2. 

We have examined this model for other values of K, 
up to K = 10, and found the same kind of transition 
in all cases. The average fraction r as a function of p 
is shown in Figure || for N = 10 5 and several values of 
K, calculated over 10 3 realizations. It is seen that the 
critical randomness p c decreases as K grows. This is due 
to the increment in the number of long-range contacts per 
element, as a consequence of the higher connectivity of 
each site. On the other hand, the value of r at p = 1, ri, 
grows with K and approaches the level expected for the 
original version of the model, r* = 0.796 . . . |1^JT6| . The 
insert in Fig. || displays the dependence of p c and of the 
difference r* — r% with K. Though this plot covers less 
that one order of magnitude in the .ftf-axis, the results 
suggest power-law decays for both quantities, p c ~ K~ p 
with p ~ 2.5 and r* — n ~ K~ a with er ~ 1.5. 



1.0 



0.8- 



0.6 



0.4- 



0.2 




0.0 



&10- 1: 5 



3 2. 



-"i — 
10 = 



10" 



10" 



10" 



In summary, our numerical analysis shows that the 
regimes where the rumor is bounded to a finite neigh- 
borhood of the initially infected site and where it affects 
a finite fraction of the population are separated by a well- 
defined transition at a finite randomness p c of the under- 
lying small-world network. The results strongly suggest 
the presence of a critical phenomenon at p c . Quite inter- 
estingly, the transition is also found if, instead of building 
a frozen small- world network as above, a dynamic small- 
world is considered ]To[j . In this case, each contact of an 
I-element is established with one of its 2K nearest neigh- 
bors with probability 1 — p, and with a randomly chosen 
site with probability p. The transition observed for these 
dynamic small-worlds is of the same type as on static 
small- world networks, but presents quantitative differ- 
ences. For K = 2, for instance, the critical randomness 
shifts to p c w 0.07. Details will be discussed elsewhere 
0- 

The geometrical properties of small-world networks 
have been shown to exhibit a cross-over from ordered 
to random behavior at p ~ -/V -1 , which implies a van- 
ishing critical value of p for asymptotically large systems 
In contrast, we have here presented evidence that a 
dynamical process occurring on such structures exhibits 
critical behavior at a finite value of the randomness. This 
difference suggests that an explanation for the origin of 
the transition in geometrical terms only is unsuitable, 
and specific dynamical properties must be taken into ac- 
count. Very recently, evidence of critical behavior at 
finite small-world randomness has been reported for a 
purely epidemiological model, whose evolution rules are 
substantially different from those considered here |jTl| . 
In particular, the epidemiological model allows for the 
transformation R — > S, giving rise to a closed disease cycle 
(SIRS) through recovery. Moreover, the transformation I 
— > R — > S is fully deterministic. The critical phenomenon 
found in that case is a transition to global synchroniza- 
tion of local disease cycles, whereas in our propagation 
process we have a kind of percolation phenomenon [||,^| . 
In spite of these basic differences, the dependence of the 
respective order parameters on the randomness p is strik- 
ingly similar, even if compared quantitatively. This sim- 
ilarity calls for further investigation in order to identify 
and characterize the whole class of small- world dynami- 
cal processes with critical behavior at finite randomness, 
and to give an analytical description of such behavior. 

Enlightening discussions with M. Kuperman and G. 
Abramson are gratefully acknowledged. 



FIG. 3. Average fraction r of refractory elements at the end 
of the evolution as a function of the randomness p for several 
values of K, and N = 10°. The insert shows the critical 
randomness p c (full dots) and the difference r* — n (empty 
dots; see text for definitions of r* and n) as a function of K. 
Straight lines in the insert correspond to power laws ~ K~ 15 
and K~ 2 - 5 . 
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